Ultra-long-range giant dipole molecules in crossed electric and magnetic fields 
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We show the existence of ultra-long-range giant dipole molecules formed by a neutral alkali ground 
state atom that is bound to the decentered electronic wave function of a giant dipole atom. The 
adiabatic potential surfaces emerging from the interaction of the ground state atom with the giant 
dipole electron posses a rich topology depending on the degree of electronic excitation. Binding 
energies and the vibrational motion in the energetically lowest surfaces are analyzed by means of 
perturbation theory and exact diagonalization techniques. The resulting molecules are truly giant 
with internuclear distances up to several jim. Finally, we demonstrate the existence of intersection 
manifolds of excited electronic states that potentially lead to a vibrational decay of the ground state 
atom dynamics. 



PACS numbers: 31.50.-x, 32.60.+i, 33.20.Tp, 33.80.Rv 
INTRODUCTION 

Since the discovery of Bose-Einstein condensates [TJ [5J 
the knowledge of ultracold atomic and molecular systems 
has grown with a breathtaking speed. Experimentally, 
one can extensively control the external motion of the 
atoms by designing and switching between almost arbi- 
trarily shaped traps [3HS] and the strength of the interac- 
tion among the atoms can be tuned by magnetic or opti- 
cal Feshbach resonances [SHE]. This has not only lead to 
a plethora of novel observations on the many-body be- 
havior of ultracold atomic ensembles [7J but also raised 
new ultracold few-body physics such as weakly bound 
diatomic molecules [S] or Efimov states of timers [TO] 
in the universal regime close to the dissociation thresh- 
old. A particularly striking new species are the weakly 
bound ultra-long-range diatomic molecules composed of 
a ground state and a Rydberg atom whose existence 
has been predicted theoretically a decade ago [TT] and 
which have been discovered experimentally only recently 
[12] . The molecular Born-Oppenheimer potential energy 
curves, which are responsible for the atomic binding, 
show for these species a very unusual oscillatory behavior 
with many local minima. The latter can be understood 
intuitively and modeled correspondingly as the interac- 
tion of a neutral ground state atom with the Rydberg 
electron of the second atom. The equilibrium distance 
for these molecular states is of the order of the size of the 
Rydberg atom and the vibrational binding energies are in 
the MHz to GHz regime for principal quantum numbers 
n « 30 — 40. Recently, the impact of magnetic fields on 
these ultra-long-range molecules [13] and the formation 
of Rydberg trimers and excited dimers by internal quan- 
tum reflection [T3] has been explored. Moreover, it it 
has been shown how the electric field of a Rydberg atom 
electron can bind a polar molecule to form a giant ultra- 
long-range stable polyatomic molecule [H] [16] and how 
Rydberg macrodimers are formed via long-range inter- 
action between two atoms excited into high-n Rydberg 



states [17] . 

An exotic species of Rydberg atoms in crossed electric 
and magnetic fields are the so-called giant dipole states 
(GDS), which have been explored theoretically [T8"rf5T] 
and experimentally [55] [53] firstly in the 1990s. Oppo- 
site to the usual Rydberg states, the GDS are of decen- 
tered character and possess a huge electric dipole mo- 
ment. More precisely, in ref.[5T] it was shown that the 
total potential of the electronic motion possesses a gauge 
invariant term which leads to an outer potential well that 
supports weakly bound decentered states. The mathe- 
matical origin of these effects is the non-separability of 
the center of mass and electronic motion in the presence 
of the external fields [55] [55] : translation symmetry and 
conservation of the total momentum in field-free space 
is replaced by a phase space translation symmetry and 
the conservation of the total pseudo momentum. Ap- 
plications to the Positronium atom have demonstrated 
that metastable matter-antimatter states with a lifetime 
of many years can be formed |26) . More recently, giant 
dipole resonances of multiply excited atoms in crossed 
fields have been shown to exist and the corresponding 
electronic configurations as well as their stability have 
been analyzed (27J [28] . 

In the present work we combine the concepts of atomic 
giant dipole states and field-free ultra-long-range di- 
atomic molecules: We show the existence of ultra-long- 
range giant dipole molecules emanating from giant dipole 
states. They exist in a variety of different configurations 
with simple to complex three-dimensional potential en- 
ergy surfaces such as Gaussian, elliptical or toroidal wells. 
For higher excited states the potential energy surfaces of 
energetically neighboring states come close in configura- 
tion space and form higher-dimensional seams of avoided 
crossings that could lead to rapid decay processes of vi- 
brational wave packets. 
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MODEL HAMILTONIAN 

We consider a highly excited hydrogen atom interact- 
ing with a ground state neutral perturber atom (we will 
focus on the 87 Rb atom here) in crossed static homoge- 
neous electric E and magnetic B fields. The correspond- 
ing Hamiltonian reads 



H = 



A. 

2m„ 



V„ 



(1) 



where 'n' labels the neutral perturber and where the first 
term is the kinetic energy of the perturber atom followed 
by the giant dipole Hamiltonian of the hydrogen atom 
and the interaction term between the GDS and the neu- 
tral perturber. The hydrogen atom in crossed external 
fields has been discussed in detail in ref. [21]. There 
it was shown that the giant dipole Hamiltonian can be 
transformed into an effective single particle problem in 
a magnetic field in the presence of a generalized poten- 
tial V(r; {E, B}) = ^(K - eB x r) 2 - ^ which para- 
metrically depends on the external fields and the pseudo 
momentum K(E, B) and contains both the motional and 
external electric field Stark terms. 
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with ii = m ° mp , o = ( m ° . " p ] e where e, to., m„, M 
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are the electron charge and mass and the proton and 
atomic mass, respectively, r, p represent the coordinate 
and canonical momentum of the Rydberg electron. For 
sufficiently strong electric fields the potential V exhibits 
an outer well containing many bound GDS which are de- 
centered from the proton at distances of about 10 5 ao, ao 
being the Bohr radius. The latter leads to a huge electric 
dipole moment typically of the order of many ten thou- 
sand Debye for strong electric and magnetic laboratory 
field strengths. Expanding V(r; {E, B}) up to second or- 
der around the minimum of the outer well and performing 
a corresponding gauge centering [201 121] we arrive at the 
second order giant dipole Hamiltonian which represents 
a charged (effective) particle in a magnetic field and an 
anisotropic 3D harmonic potential, 
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where the frequencies u> x = yj |(|j7 + ^), u y = 

\J)l(£m ~ Jf) and Wz = l/VlM^ol characterize the 
anisotropy of the outer well. In this representation, r, p 
denote the electronic variables with respect to the outer 
minimum (xo,0,0) with xq ~ —[K\/B. Our working 
Hamiltonian therefore reads 



H = 



Pi 

2to„ 



— (p-^Bxr) 2 
2// P 2 ' 



+ ^ + (Vy 2 + |Vz 2 + K,oo(r,r n ), (4) 



where K.gd again represents the interaction of the neu- 
tral perturber atom with the giant dipole electronic Ryd- 
berg state. For deeply bound states in the outer well the 
electron possesses a low kinetic energy and it is legitimate 
to describe the interaction with the neutral perturber by 
a Fermi-type pseudopotential jTTJ [2U] , namely, a s-wave 
contact potential 



K.QD(r,r ) = 2 7 r^ T [£ ; (r)]<^(r-r„). 



(5) 



Here A T (k) is the energy-dependent triplet s-wave scat- 
tering length for electron collisions with the ground state 
Rb atom. r n denotes the position of the neutral per- 
turber with respect to the minimum of the outer well. 
The electron wave number k is provided by the kinetic 
energy of the Rydberg electron when it collides with the 
neutral perturber [IT] . 



METHODOLOGY 

In order to solve the eigenvalue problem associated 
with Hamiltonian Q we adopt an adiabatic ansatz for 
the neutral ground state and the giant dipole atom. We 
write the total wave function as ^(r, r„) = cf>(r„)ip(r; r„) 
and yield 

[H GD + K, GD (r, r n )]V>»(r; r n ) = e,-(rJV>i(r; O, ( 6 ) 

r^2 



2to„ 



e*(r„) 



0UO = *M(r n ), 



(7) 



where tp describes the electronic wave function of the de- 
centered GDS in the presence of the neutral perturber 
for a given position r„ and <f> determines the vibrational 
state of the neutral perturber. To calculate the poten- 
tial energy surfaces (PES) Ci(r n ) we expand tpi(r;r n ) in 
the eigenbasis of H GD , i.e., ^(rjrj = £j ; Cj(r n )xj( r ) 
with H GU Xj{r) = SjXj(j), and solve the corresponding 
eigenvalue problem associated with eq. ^ using stan- 
dard numerical techniques for the diagonalization of Her- 
mitian matrices. In |21j it was shown that the eigenen- 
ergies and eigenfunctions of H Gr> are determined by the 
three quantum numbers n_, n +1 n z = 0, 1, 2, ...[xj(r) = 
Xn_n + n z (r)], where e„_, n+ ,„ z = huj-(n- + \)+huj + (n + + 



) + (n z 



with u>± = \u>z 



l 

\/2 l 



■ oj 2 ± sgn(w 2 



l 

2 

w y 2 )^(w 2 + w 2 + w 2 ) 2 - ^w 2 ] 1 / 2 and uj c = -qB/fi. To 
ensure convergence for our numerical approach we vary 
the number of orbitals associated with the quantum num- 
bers n_ , n + , n z independently finally achieving a relative 
accuracy of 10 -5 for the energy. To do so for the energet- 
ically lowest 15 excitations a basis set of approximately 
1500 states is needed. In addition to the numerically 
exact treatment we determine the PES in first order per- 
turbation theory, leading to 



ef(rJ= £j +2^ T [fc(rJ]|x,(r n )| 2 



(8) 
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From eqs. Q-Q we can already deduce some sym- 
metry properties of the states "J, <f> and the ener- 
gies e(r). If P r ,r n denotes the parity operator that 
transforms (r, r n ) — » (— r, — r a ) we have [H, P T ,r n ] — 
[K,GD(r, r n ), P r ,r n ) = 0. This means that the states ty, ijj 
and c/> are parity (anti) symmetric and the PES are sym- 
metric, i.e., e(±r„) = e(r n ). 

Throughout this work we use the exemplary field con- 
figuration E = 50 V/cm and B = 2.35 T. This gives 
giant dipole level spacings of w_ = 223 MHz, oj + = 
413 GHz and uj z — 1.35 GHz. For such fields, it is jus- 
tified to neglect all interaction terms due to the induced 
dipole of the Rb atom. 



MOLECULAR GROUND STATE POTENTIAL 
SURFACES 

In first order perturbation theory the molecular giant 
dipole ground state PES is given by 



c 000 



(r„) =£ 000 + 27rA T [/c(r n )]Coe 



2 i 2 2 

—ax„—oy n —cz- Q 



(9) 



where the constants Co, a, b, c > are given in 21j. This 
potential represents a single well 3D potential and since 
Ar[fc(r n )] < Vr„ the minimum is located at r n min = 
with a depth of Ae Q00 /h = [e 000 - e£ 00 (r nmin )]//i = 591 
MHz. In the present case we have A T [fc(r n )] sa A T (x^ + 
y^,z^) and a rs 6, which means that the giant dipole 
ground state PES possesses an approximate </? n -rotational 
symmetry and (L z )„ is an approximately conserved quan- 
tity. A harmonic fit V h ( r n) = ^ m a UJ x x l + \ m ^ uo yyt + 
^m D 0J^z^ around this minimum provides harmonic oscil- 
lator states with a level spacing of u) x = u> y = 39.2 MHz 
and uj z = 12A MHz. The numerical exact PES is pre- 
sented in Fig. [TJa). The employed basis set consists of 
N- =30, N+ = 1, N z = 50 states. As in perturba- 
tion theory the potential possess a minimum located at 
r n.min = and is iy9 n -rotationally symmetric. A harmonic 
fit provides here a depth of Ae oo/^- = 788 MHz and level 
spacing of u> x = io y = 20 MHz, us z — 6.5 MHz. 

In Fig. [ljb) we show a comparison between the per- 
turbative and the exact ground state potential for y n = 
z n = 0. In addition, the harmonic fitting curves and 
the parameters uj* and w° are presented as well. Com- 
paring the analytic solution eq. (|9| with the numerically 
exact result, one finds that for the exact PES the depth 
possesses a 165 MHz larger value while the width of the 
well increases by a factor of two. Furthermore, for the 
perturbative PES the harmonic approximation is valid 
in a region between x„ ~ ±200ao, while for the exact 
potential it is valid up to x n rs ±800ao- Hence, while 
providing a good qualitative prediction, the perturbative 
approach cannot be used to discuss quantitative details. 
The deviation between the exact and the perturbative 
result can be understood by analyzing the coupling of 
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FIG. 1. (a) Ground state potential energy surfaces calculated 
via exact diagonalization. The employed basis set consists 
of iV_ = 30, iV+ = 1, 7V Z = 50 states. The energy scale is 
given in GHz. (b) Comparison between perturbative (dashed- 
dotted line) and exact (solid line) ground state potential for 
y„ = z n = 0. In addition, the harmonic fitting curves (thin 
lines) together with the corresponding trap frequencies 
and wj x are shown, (color online) 



the giant dipole groundstate to excited states of the un- 
perturbed system. For first order perturbation theory to 
hold, these couplings need to be much smaller than the 
energetic separation of the involved levels. In our case, 
this amounts to the requirement 



7n_ 



:= 4 



(000|K,G D (r,r n )|72_n + n z ; 
(r n ) - e%_ n+n „(r D ) 



c ooo 



< 1 (10) 



where n_ , n + and rt z label the quantum number of the 
excited state. In our case, we find, e.g., 7100 ~ 4.5 and 
7002 ~ 0.1. Hence, it is not surprising that first order 
perturbation theory yields qualitatively but not quanti- 
tatively reliable results. This is in contrast to the trilo- 



bite systems [TT] where degenerate first order perturba- 
tion theory within a given Rydberg n-manifold provides 
satisfactory results (depending on n, the splitting of ad- 
jacent manifolds is in the GHz-THz regime compared to 
cj_ = 223 MHz in our case). 

We remark that the vibrational states obtained by 
the potential given in Fig. 1 are localized at internu- 
clear distances in the range of 10 5 ao- To our knowledge, 
these molecules belong, together with recently investi- 
gated Rydberg macrodimers [17] , to the largest diatomic 
molecules ever predicted. 



POTENTIAL SURFACES OF EXCITED STATES 

In perturbation theory the PES for the lowest excita- 
tions are determined by the states |n_00), n_ = 1,2,3. 
Introducing cylindrical coordinates p n , ip a , z n and ex- 
tracting the dominant terms in perturbation theory, the 
PES are approximately given by 



C 7l_00 



(r n ) « e n _ 00 + 27rA T [k(v a )}C n ,e^- c ^pl n in) 



Due to the weak dependence of A T [k(r n )] on k in our case, 
this well represents a t/? n -rotationally symmetric torus 



with minima at r 
and depths Ae ra 
given by 

>-) 
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P mi „ ' « y/n-/a, 

Ae n _o *2K\A T [k(r£d)}\^(—) n ~ 

e n - \ a / 



(12) 



In Fig. [2]ja) the exact PES for the first excited state 
is shown. The toroidal structure as predicted by per- 
turbation theory is clearly visible. Fig. [2] (b) shows a 
2D intersection for z n = of the same PES: A "Mex- 
ican hat" like potential well is obtained with the one- 
dimensional manifold of the minimum lying on a cir- 
cle. For sufficiently small displacements from the min- 
ima the PES can be described by a harmonic approxi- 
mation, i.e., a decentered oscillator potential V h (r D ) = 
|m n Wp(p n — p min ) 2 + \m I jJ\z^. One therefore arrives 
at a free rotational motion of the Rb-atom in the in- 
direction. In perturbation theory the frequencies w p and 
lo z are given by 



'27ro|A[A;(ri^)]|C, 



,(n-) 



1 nc\A[k(r£-l)]\C n _ 



(14) 



The corresponding eigenfunctions are given by the prod- 
uct A(p n )e imVa 4>i(z tL ), m G Z and I e N , where fa de- 
notes the l-th eigenfunction of the harmonic oscillator 
and A(p n ) are, apart from a Gaussian, the biconfluent 
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FIG. 2. First excited potential surface. To achieve conver- 
gence we use a set of iV_ = 10, N+ = 1, N* = 30 basis 
functions, (a) 3D representation and (b) 2D intersection for 
z n = 0. The energy scale is given in GHz. (color online) 



Heun functions [30] ■ In table [I] the parameters p mln , w p 
and w a for the first three excited PES, both exact and 
perturbative, are listed. We see that the largest differ- 
ence between perturbation theory and the numerically 
exact result is obtained for the p min and ui z parameters 
with a deviation of around 50% for the first excitation 
and 20% for the third excitation. However, for uj p and 
Ae ra _oo both results are more comparable with a maxi- 
mum deviation of 20%. For all quantities the exact and 
perturbative results are more comparable with increasing 
excitation. 
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exc. 1 


exc. 2 


exc. 3 




pt 


ex 


pt 


ex 


pt 


ex 


Pmin / «0 


450 


825 


637 


971 


785 


1050 


Ae„_ 00 / MHz 


221 


170 


162 


155 


134 


127 


Up 1 MHz 


33 


27 


29 


25 


26 


24 


uj z 1 MHz 


1.9 


0.86 


1.7 


0.89 


1.5 


1.2 



TABLE I. Parameters p min , Ae n _oo, uj p and oj z for the first 
three excited potential surfaces. For the numerically exact 
potentials (ex) the parameters are taken from a harmonic fit- 
ting. For the perturbative potentials (pt), the parameters are 
extracted from eqs. ( 12 1-( 14 1 . 



AVOIDED CROSSINGS OF POTENTIAL 
SURFACES 
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While for the first few excitations first order pertur- 
bation theory provides a reasonable prediction of the 
qualitative behavior of the PES, for higher excitations 
first order perturbation theory is not capable of describ- 
ing even the qualitative behavior of the PES. For ex- 
ample, in Fig. [3] we show the PES for the fifth excita- 
tion. In addition to the toroidal well described above 
we get two new 3D elliptical potential wells centered at 
x n = y n — 0j z n ~ ±6500ao. As a consequence the former 
rotationally symmetric global minimum represent now a 
local minimum and the two new global minima are the 
centers of these elliptical wells. The bound states in the 
toroidal well become metastable and can in principle de- 
cay into bound vibrational states in the elliptical wells. 

For the fifth and sixth excited PES a harmonic ap- 
proximation around the global minima yields a vibra- 
tional level spacing of approximately 20 MHz. The 
additional elliptical wells arise due to avoided crossing 
with neighboring PES. To confirm this, Fig. [4^a,) shows 
intersections for the fifth, sixth and seventh PES for 
x n = 0, y n = 200eto. Avoided crossings of the PES are en- 
countered, e.g., at z w ±8.5 x 10 3 ao and z ~ ±5.5 x 10 3 ao 
for the fifth and sixth excited PES. 

A more global view of the geometry of the higher ex- 
cited PES and their avoided crossings is provided with 
Fig. Kb). It shows the PES obtained via perturba- 
tion theory, i.e., eg* 00 (r o ), eg 00 (r n ) and e% 01 (r n ) for fixed 
x n = 0. The potential eooi(r„) is obtained by raising the 
quantum number n z and is given by 

eg 01 (r n ) = eooi + 2TTA T [k(r n )}Cz 2 e- a p'- cz ". (15) 

In contrast to £g 00 (r n ) and £g 00 (r n ) this potential curve 
possess two minima at z = ±5542ao. Because the 
depth of these minima (« 450 MHz) is larger than 
the level spacing of the unperturbed giant dipole levels 
[(e oi - £ 6 oo)A « 9 MHz, (e 001 - e 50 o)/fc ~ 250 MHz] 
the PES egQ 1 (r n ) intersects neighboring PES. Due to the 
coupling between the different PES this leads to avoided 



FIG. 3. Potential surface for the fifth excitation. The addi- 
tional elliptical wells arise due to avoided crossings with ad- 
jacent potential surfaces. The energy scale is given in GHz. 
(color online) 

crossings for the exact PES and consequently novel ge- 
ometries of the potentials. In the vicinity of the avoided 
crossing the adiabatic approximation fails and we expect 
a strong rovibronic interaction mixing different electronic 
giant dipole states. As a consequence fast decay processes 
of wave packets probing the seam of the avoided cross- 
ings will take place. Finally we note that the radiative 
lifetimes of the excited states are for our chosen parame- 
ter values of the order of several days and therefore much 
longer than the typical vibrational frequencies in the ex- 
cited PES: it should therefore be possible to probe the 
vibrational dynamics belonging to the complex geometry 
of the molecular PES of the excited giant dipole states. 

We remark that the occurence of avoided crossing will 
be even more prominent in heavier giant dipole systems. 
Considering rubidium instead of hydrogen, for example, 
reduces the smallest level spacing by m p /M Rb s» 10~ 2 , 
which yields cj_ =2.5 MHz. The depths of the potential 
surfaces, on the other hand, do not change significantly. 
As a result, already the ground state shows avoided cross- 
ings with neighboring potential surfaces. 

CONCLUSION 

Bringing together the concepts underlying atomic gi- 
ant dipole states in external fields and ultra-long-range 
molecules, we demonstrated the existence of ultra-long- 
range giant dipole diatomic molecules. In particular, the 
exotic atomic state underlying these molecules give rise 
to novel properties such as a plethora of different quan- 
tum states with complex three-dimensional energy land- 
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FIG. 4. (a) Intersections of the fifth, sixth and seventh excited 
potential surfaces for x n = 0, j/„ = 200ao are shown. Avoided 
crossings of the potentials are encountered. The spacing be- 
tween the fifth and sixth surface at the avoided crossings is 
approximately 9 MHz. Subfigure (b) shows 2D intersections 
of the perturbatively calculated PES for x a = 0. (color online) 



scapes and rich rovibrational dynamics. The resulting 
molecules possess very large rovibrational bound states 
at internuclear distances in the range of 10 5 ao. 

For their experimental preparation the 'best of two 
worlds' has to be combined. The preparation of giant 
dipole states is known to be possible starting from 'tra- 
ditional' Rydberg states in magnetic fields and applying 
a sequence of electric field switches which brings the elec- 
tron into low- lying outer well states [31] ■ Driven radio- 
frequency transitions in the outer well as an additional 
tool might help to prepare definite outer well states. 
Starting from these, one could overlap the GDS with a 
dense cloud of ultracold rubidium atoms and use radio- 
/ microwave induced transitions to form the envisaged 
giant dipole molecular states. One of the main differ- 
ences to standard cold atom experiments is certainly the 



regime of field strengths necessary to address the giant 
dipole states, which corresponds to strong static mag- 
netic and electric fields. 
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